Classical Strongly Coupled QGP: 
VI. Structure Factors 

Sungtae ChcQ and Ismail Zahecfl 

Department of Physics and Astronomy 
State University of New York, Stony Brook, NY 11794-3800 

Abstract 

We show that the classical and strongly coupled QGP (cQGP) is characterized by a multiple of 
structure factors that obey generalized Orstein-Zernicke equations. We use the canonical partition 
function and its associated density functional to derive analytical equations for the density and 
charge monopole structure factors for arbitrary values of V = V/K, the ratio of the mean potential 
to Coulomb energy. The results are compared with SU(2) molecular dynamics simulations. 
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I. INTRODUCTION 



High temperature QCD is expected to asymptote a weakly coupled Coulomb plasma 
albeit with still strong infrared divergences. The latters cause its magnetic sector to be non- 
perturbative at all temperatures. At intermediate temperatures of relevance to heavy-ion 
collider experiments, the electric sector is believed to be strongly coupled. 

Recently, Shuryak and Zahed |l| have suggested that certain aspects of the quark-gluon 
plasma in range of temperatures (1 — 3) T c can be understood by a stronger Coulomb in- 
teraction causing persistent correlations in singlet and colored channels. As a result the 
quark and gluon plasma is more a liquid than a gas at intermediate temperatures. A liquid 
plasma should exhibit shorter mean-free paths and stronger color dissipation, both of which 
are supported by the current experiments at RHIC ^]. 

To help understand transport and dissipation in the strongly coupled quark gluon plasma, 
a classical model of the colored plasma was suggested in [3j. The model consists of massive 
quarks and gluons interacting via classical colored Coulomb interactions. The color is as- 
sumed classical with all equations of motion following from Poisson brackets. For the SU(2) 
version both molecular dynamics simulations 3j] and bulk thermodynamics were recently 
presented 4] including simulations of the energy loss of heavy quarks Q]. 

In this paper we follow up on our recent equilibrium analysis of the bulk thermodynam- 
ics (4] to the static structure factors. In section 2 we define the energy functional for the 
cQGP. In section 3 we derive generalized Ornstein-Zernicke equations for the pair corre- 
lation functions. In section 4 we show that the cQGP supports multiple structure factors 
that measure a variety of colored correlations. Each structure factor obeys a generalized 
Ornstein-Zernicke equation. In section 5 we introduce the Debye-Huckel-hole potential for 
the cQGP. In section 6, we use Debye charging process to derive analytical expressions for 
the lowest two structure factors in the cQGP for arbitrary T. In section 7, we construct 
numerically the lowest two structure factors using molecular dynamics simulations and com- 
pare them with our analytical results for values of T in the liquid phase. Our conclusions 
are in section 8. Appendix A is added to streamline our conventions for the SU(2) color 
charges. 
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II. FREE ENERGY FUNCTIONAL 



We consider the canonical partition function of a single species, either quarks and gluons, 
at finite temperature 1/(3 = T and in the presence of an external scalar source ip 

ZnM = M / II ex P O 3 / drdQn(r, Q)^(r, Q)) 

x exp ( - ^ f drdr'dQdQ'n(r, Q) Q ' Q n(r', Q')) (II.l) 
V 2 An J \r — r'\ / 

The color charges are treated classically and we refer to [^SIs] for further details regarding 
the nature of the measure. Here, we have defined 



N 

n' 



r,Q) = '£ i 6(r-r i )6(Q-Q i ) (11.2) 

i 

The generalization to many species is straightforward. The associated Coulomb parameter 
is 



9 2 PQ 



2 (II.3) 



47r a W s 

where C2 is the quadratic Casimir (= £V Q 2 /(N 2 — 1)) and aws is the Wigner-Seitz radius 
^ na ws/3 — 1/n. For small T, Eq. fill, lj) behaves as a screened but weakly coupled gas, 
while for intermediate values of T, Eq. (III. ID describes a liquid 0, (|. At large values of T 
Eq. fill. 1 p yields a solid as a ground state. From now on, the canonical charge of g 2 /4-7T will 
be set to 1 for simplicity, and will be restored in the final parameters by inspection. 

The static correlations both in space and in phase space associated with Eq. flll.lj) are 



involved and will be the subject of most of this paper. For that, we note that Eq. fill. 1 j) 
yields the free energy generating functional 



?nW\ drdQn ( - 1 \r,Q)(^\n(X 3 n^(r 1 Q)) - 1^ 

- J dr^r, Q)n {1 \r, Q) + f c (n {1 \r, Q)) 

+ - [ drr'dQdQ'nW (r, Q) Q Q , (V , Q') (II.4) 
2 J \r — r I 



Here we have set 
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N 

(1) 



r, Q) = (n(r, Q)) = (r - r,)5(Q - QJ) (II.5) 



as the expectation value with the averaging carried using Eq. (III. II) . The second contribution 
in Eq. (111.41) is the ideal classical contribution following from the measure in Eq. (111.11) using 
the asymptotic Stirling formulae. The third contribution is the excess free energy functional. 
J- c is the connected free energy that sums up the second and higher cumulants of n(r, Q) 
from Eq. (111.11) . We note that for zero scalar source '0 = 0, 

F N [0] =Fid + Fex (116) 

where the first contribution is the classical ideal part and the second contribution the excess 
part. 

III. ORNSTEIN-ZERNICKE EQUATIONS 



To quantify the static interactions between pairs of particles in (IH.lj) we define 
1 5 2 F N 



= (n (1) (r, Q)n {1 \r', Q')h{r -r',Q Q') + n (1) (r, Q)5(r - r')5{Q - Q'fj 

(III.l) 

with h the pair correlation function for — 0. The pair correlation function is invariant 
under space translation and color rotation. Generically 

h(r -r',Q-Q') = ^(J2 s ( r ~ r ^ r ' ~ r MQ - QMQ' - Qj)) (ni-2) 

i¥=3 

The direct correlation function cp follows from the excess free energy (111.61) through 

--c D (r -r',Q Q>) = = + j (111.3) 
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which plays the role of a correlated potential, will be used below as a renormalized 
Coulomb potential in the liquid phase. It also obeys the identity 



= -m = ji O " r<) w " Q,) ~ c ° [r ~ r ' Q ' Q '» l (IIL4) 

Using the chain rule, 



dr"dQ" 6 ^ Q) ^V^j - tf r _ r 'U(Q - Q') (III 5) 

we obtain 



h(r-r',Q Q') = c D (r-r',Q-Q')+n J dr"dQ"h(r-r" , Q-Q")c D (r"-r', Q' Q") (III.6) 

which is the Orstein-Zernicke equation that ties the pair correlation h to the direct correla- 
tion or the pair potential cp. For a uniform plasma ( 1III.6j) unfolds agebraically in momentum 
and color space using 



h(r -r',Q qf) = f dke ik < r ~^ ]T 2 I±l\ ll {k)P l {Q ■ Q') 
J i 71 

c D (r -r',Q Q') = J dke ik < r '^ ]T 2 l±l CDl (k)P l (Q ■ Q') (III.7) 



Thus 



hj(fc) = c Dl (k) + nhi(fc) c Di (fc) (III.8) 

which holds for each partial waves. (1III.8P are the generalized Orstein-Zernicke equations 
for each color partial wave of the SU(2) colored Coulomb plasma. 

IV. STATIC STRUCTURE FACTORS 

The statistical aspects of the colored charged particles are best captured by correlations 
in the phase space distributions. The static structure factor is defined as 



S (r - rW, Q Q ) = (5f(rpQ)5f(r'p'Q>)) 
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(IV.l) 



with formally 



f(rpQ) = S(r - Xi)S(p - Pt )5(Q - Qt) (IV.2) 

i 

and 5f — f — (f). The averaging in (llV.lj) is carried using the canonical partition function 
(IH.lj) . Color and translational invariance imply 

(f(rpQ)) = nf {p) = n(-^) 3/2 e-^ 2 / 2m (IV.3) 

which is the Maxwellian distribution for massive constituent quarks or gluons. It is readily 
shown that 



S (r-r',pp',Q-Q') = nf (p)5(r-r / )5(p-p / )5(Q-Q') + n 2 f (p)f (p')h(r-r',Q-Q') 

(IV.4) 

The reduced static structure factor 

S (fc, Q-Q') = ^ J dpdp' J dke ik < r - r,) S (r - r',pp', Q ■ Q') (IV.5) 
ties with the pair correlation function 2 j) through 

S (fe, Q Q') = 5(Q - Q') + nh(k, Q ■ Q') (IV.6) 
Its Legendre transform in the color charge reads 

So^fe) = 1 + nh(k) (IV.7) 

So the knowledge of the partial- wave structure factor <Sj(A;) yields the pair correlation c Dl 
through ffTTTT7|) and (jTVT7l) 

1 = Soi(k)- 1 + nc Dl (k) (IV.8) 
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We note that in configuration space the 1th partial wave of the static structure factor is 



Soj(r-r') = i J dp dp' J dQ P t (Q ■ Q') S (r - r', pp', Q ■ Q') (IV.9) 



Using (1III.2j) and (1IV.4j) and enforcing space translational and color rotational invariance in 



the averaging process yields 



Mr) = 6{r) + i( 5< < r ~ r ^ P tQi ■ Qi)) ( IV - 10 ) 
In particular, the two lowest static structure factors are the density structure factor 



Soo(r) = 5(r) + 6 ( r ~ r *i)) = S ( r ) + nh o(r) (IV. 11) 

and the charge structure factor 



Soi(r) = 6(r) + ^( E *( r ~ r ^ Qi ' Qi) ( IV - 12 ) 

Higher structure factors are given by fllV.lOj) as they measure the various color correlation 
content of the SU(2) strongly coupled QGP. Below, we propose both an analytical and 
numerical derivation of the two lowest structure factors flIV.111) and (jIV.12j) . 

V. DEBYE-HUCKEL-HOLE POTENTIAL 

To derive the static structure factors we will use the Debye charging procedure for a 
fixed color charge. For that, we need the Poisson-Boltzman equation for the 1-species SU(2) 
colored plasma in the presence of a colored test charge q [3| 



V 2 0(r, r', q) = -4vr ^q5(r - r') + J dQ'Q'n(r, q ) e -PQ'i<Hr,r>, q )-*(r,<i)) j 



(V.l) 



with the fixed external density profile 



n ( r i q) — n + n (^o + Ai ■ q) cos(fc • r) (V.2) 
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We note that (jVT2]) is a scalar under rigid and orthogonal color rotations HQ if the external 
parameters A; transform unitarily as _D(R)A with -D(R) the Wigner rotation in the adjoint 
representation. This fixed density causes an imposed potential 



V 2 $(r,q) = -4vr q(n(r,q) - n) (V.3) 



which is used to normalize the Poisson-Boltzman equation in ( IV. II) . We solve ( IV. II) in the 
linear approximation. For that we define the shifted potential 5(f) = (/> — $, 

(V 2 - K 2 D n(r, q) /nj 6(j)(r, r', q) « -47rq (5{r - r') - (n(r, q) - n) j (V.4) 

with n 2 D = 4ir(3nC2 the squared Debye constant. (1V.4I) is the linearized Poisson-Boltzman or 
Debye-Huckel equation for the classical colored plasma. At short separations (IV.4j) is known 



to misrepresent the hole caused by the strong Coulomb correlations. To fix that we use the 
Debye-Huckel plus hole approximation 1 1 Of] 

(V - K 2 D n(r, q)/nO) 5<j>(r, r', q) w -4vrg (5(r - r') - (n(r, q)-n))(l- 0) (V.5) 

with = 6(\r — r'\ — a) the spherical hole insertion of radius a. The mean-induced potential 
is 

*(r, q) = lim (6(f>(r, r', q) - , Q ,. ) (V.6) 

VI. DEBYE CHARGING PROCESS 

To assess the static structure factors for the classical and strongly coupled colored plasma, 
we note that the excess free energy (111.61) can be readily rewritten in terms of the pair 
correlations 

-{3F ex = I J drdr'dQdQ'n^\r, Q)c D {r -r',Q Q')n^\r', Q') (VI. 1) 



8 



In Fourier (space) and Legendre (color) space (IVI.lj) reads 



-P? ex = \f dkdQdQ'J2n2(k 7 Q)c Dl (k)nl 1 2(-k,Q') (VI.2) 

Lm 



with 



n£(fc, Q) = J dre- ik - r Yr(Q)nU(r, Q) (VI.3) 

and Y" z m a spherical harmonic for an SU(2) colored plasma. Using the partial wave form of 
the Orstein-Zernicke equations (1IV.8I) . (1VI.2[) becomes 



= ~L\ dkd Q d Q'J2 n ^Q)( 1 ~ S f (fc))n«(-fe,Q') (VI.4) 

l,m 

This shows that the quadratic change in the excess free energy caused by an external density 
profile ni m (k) is directly proportional to the Ith partial wave of the inverse of the static 
structure factor. 

The external density profile (1V.2|) changes the color Coulomb potential locally through 



flV.4p . thereby affecting the free energy. To assess the change in the latter we use the Debye 
charging procedure For that, we note that by dialing flV.2|) the free energy shifts. The 
shift can be decomposed into three parts, 

•F -^"ideal -^"imposed -^"induced 

The shift in the ideal part is set by the first term in ( 111.41) after inserting ( IV. 2 1) . The imposed 
contribution is 

imposed = J drdQ(n(r,q) - n)q ■ J d\$(r,\q) (VI.6) 
and follows from the imposed charge. Specifically, 



^ n 2 9/„.9 .,,,[,, ,cos(fc • r) cos(k ■ r') 

imposed = y? 2 (3A2 + |Ax| 2 ) / drdr' 1 (VI.7) 
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The induced free energy is 



^induced = J drdQn(r, q)q ■ j dX^(r, Xq) 
and follows from the induced but shifted screening; potential (TV6l) . 



(VI.8) 




FIG. 1: Soi(g) for T = 4,8, 12, 16 summed up to / = 1 (a) and I = 2 (b). See text. 



The integrand can be obtained by solving ( IV. 51) for 5<p(r, r', q) with the help of the Green 
function, 



V 2 - K 2 D B(\r-r' 



(J 



G(r", r-r') = -An5\r" - (r - r')) 



(VI.9) 



The method has been developed in 
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141 ] for the one-component plasma and readily 
extends to our colored plasma. For that we evaluate the reduced free energy / = (3 J- '/V 
to quadratic order in A 2 . By comparing the terms with (IVI.4I) . we can extract Sqq 1 as the 
coefficient of Aq and as the coefficient of A 2 . We find that both static structure factors 
are finite and identical in this approximation, 
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Soo^fe) = S m \k) = 1 + 2(^) 2 f dX— cos (\ — (w x - 1)) 

k J w x Xk d 



, dA— sin -— (u; A -l -2 / — — T+(w x -l,-— , 



+E(a+i)(*T f dX iWx 7 ' r ,^ -(l> " DW"* " Lit-) 

^ \ k J w x gi +1 {w x - I) X k Xk d 

- 1, \^)X-{ Wx - l, I A) _ 27»(w x - 1, ~))) (VI.10) 

A Kf) A Kd A K]j / J 

The three integral contributions are 

I r( z ^y)= dww- l g l (-w)ji(yw) 
Jo 

POO 

J X 

POO 

l+(x,y) = dzz- l gi (z)Uyz)e^~^ (VI.11) 



with ^(z) = e z z l+1 ki(z). Here ji(yz) is a spherical Bessel function, and fcj(z) a modified 
spherical Bessel function. The parameter w\ is defined as 



i 

Wx = ^l + X 3 (3T)iy (VI.12) 

The static structure factors in (IVI.lOp involve summations over multiple partial waves. The 
sums are rapidly converging as we show in Fig. 1 with / = 1 retained (a) and I = 2 retained 
(b). Here q = ka-^s is a dimensionless wave- vector. To assess the accuracy of the analytical 
method developed above for the static correlation functions in the colored Coulomb plasma, 
we now carry numerical simulations for the same structure factors using molecular dynamics 
simulations. 



VII. STRUCTURE FACTORS FROM MOLECULAR DYNAMICS 

For an SU(2) plasma, the details of the molecular dynamics simulations can be found 
in Color motion is treated as a point coordinate on a 3-sphere with a fixed radius that 
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is equal to the quadratic Casimir for SU(2). Classical stability of the colored Coulomb gas 
at short distances is achieved by using a scalar core potential of the type 



V„ 



1 1 



n r,- 



(VII. 1) 



with n = 9. The two-body interparticle colored potential is [3j, Il6| . 



V(r,Q.Q') = 9 ^[l^) 9 + Q.Q'^)] (VII. 2) 

with A setting the unit of length scale. At close packing the density is n cp = 1/A 3 . We 
choose the unit of length A so that n cp = 1. The unit of time is set by the inverse plasma 
requency r = u' 1 . In these units, the strength of the colored Coulomb potential is j^-^s 

I 

We have adopted the Verlet algorithm in [l5j to integrate the equations of motion for 
a system composed of 108 particles. The particles are confined in a box and surrounded 
by images via periodic boundary conditions. The simulations are carried in a fixed volume 
l/(n\ 3 ) = V/(NX 3 ) = 2.72 as in Q. The Wigner-Seitz radius a ws = (Ann/3)- 1 / 3 is 0.866A. 
With these parameters, the interparticle interaction strength is set by the Coulomb constant 

r. 

We first measure the particle radial distribution function g(r) = h (r) as a function of 
T. g(r) measures the probability of finding two particles between r and r + Ar, 



g(r) = ho(r) = ^( J>(r - r y )) (VII. 3) 

In Fig. [2] we show g(r) versus r for different Coulomb couplings Y = 2.2, 6.6, 12.8. The larger 
T the larger the size of the Coulomb hole surrounding each colored Coulomb particle. Also, 
the larger T, the higher the peak, the tighter the Coulomb packing. 

The radial distribution function in Fig. [2] for the SU(2) colored Coulomb plasma appears 
overall similar to the one observed for the one component plasma (OCP). Although our 
colored particles attract for color antiparallell charges, they overall statistically repel due to 
the larger color repulsive orientations. The difference with the OCP is best seen by taking 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 
r 

(b) 





FIG. 2: Radial distribution function for T =2. 2(b) , 6.6(c) and 12.8(d). See text. 



the Fourier transform of (1VH. 3j) which is the / = density structure factor Soo(<?) 



Soo(fc) = ^(|n fc | 2 ) 



(VII. 4) 



with 



n fc 



N 



(VII. 5) 



i=l 



Fig. [3] shows the behavior of S o(<?) versus the dimensionless wave-vector q = awsk- The 
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nonvanishing of Soo at the origin reflects on the coupling to the sound mode. In the static 
and long-wavelength approximation it is just 



Soo(fc) 



k 2 

c 2 s k 



(VII. 6) 



with c| = (dP/dpjx the isothermal squared speed of sound, with P the pressure and p the 
mass density. Since k is a multiple of 2ir/L because of the finite cubic box L x L x L, only 
about a dozen points were accessible numerically. Since L w N 1 ^ 3 , we need to increase the 
number of particles in the box to smoothen out the structure factor in momentum space. 




FIG. 3: Static Structure factor Soo(g) versus q for T = 2.2,6.6, 12.8. See text. 
In Fig. HI we show the I = 1 or charge structure factor 



Soi(k) = l(\ Pk \ 2 ) (VII. 7) 



N 

with 



N 



Unlike S o which correlates a pair of scalar densities, S i correlates a pair of charge densities. 
In the OCP plasma both correlators are identical. They are not in the SU(2) colored 
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Coulomb plasma. In the long wavelength approximation, the static density structure factor 
is saturated by the plasmon mode 




FIG. 4: The static structure factor Soi(<z) for T = 2.2, 6.6, 12.8. See text. 

Our analytical result for Soi(fc) in (1VI.10I) is in agreement with the molecular dynamics 
simulations for the charged correlator (I VII. 71) . Our analytical result for Soo(^) is identical 
with Soi(k). It differs from the molecular dynamics simulation results for small momenta 
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since the sound mode drops out of the Debye-Huckel colored potential on which our charging 
process was based. The contribution of the sound mode is additive at small momenta, and 
drops out at large momentum due to damping through the shear viscosity. 

VIII. CONCLUSIONS 

The strongly coupled SU(2) QGP is characterized by a number of static correlators in 
phase space with color treated as a classical 3-vector on S 3 with a radius fixed by the 
second Casimir. Space translational invariance and color rotational invariance yields multiple 
structure factors characterizing color correlations with color charges sourced by Legendre 
polynomials. Each structure factor obeys a generalized Ornstein-Zernicke equation. 

To evaluate analytically these multiple structure factors, we have made use of the Debye 
charging process and the linearized Poisson-Boltzman equation in line with linear response 
theory. We have derived explicit relations for the two lowest structure factors, ie I = 0, 1 
which corresponds to the density and charge structure factors. 

To check the validity of the linear response analysis, we have numerically extracted the 
density and charge static structure factors using SU(2) molecular dynamics simulations. 
Modulo the sound mode, both analytical structure factors compare favorably with the nu- 
merical results. The current analysis extends to higher multipoles, ie / = 2, 3, ... and gener- 
alizes to higher color SU(N > 2) groups. 

The static structure factors play an important role in characterizing the correlations in 
the colored SU(2) QGP at intermediate and large values of the coupling coupling T. They 
also enter in the assessment of transport parameters at strong coupling. The results will be 
presented elsewhere. 

The current classical and strongly coupled SU(2) colored Coulomb plasma can be ex- 
tended to several species to account for gluons, quarks and antiquarks [3]. The effects of 
quantum mechanics being a renormalization of the constituent parameters such as the mass 
and charge. It will be interesting to see whether a quantum phase space formulation of 
QCD is achievable through the background field formulation in a way that allows for the 
introduction of colored static structure factors. 
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APPENDIX A: SU(2) COLOR CHARGES 

The explicit representation of the classical color charges is js, 9] 



Q 1 = cos 0i J J 2 - tt 2 , Q 2 = sin0 n / J 2 - tt 2 , Q 3 = tti (A.l) 



with J 2 the quadratic Casimir q 2 = J2 a c Q a Q a - The measure in the SU(2) phase space 
can be set to 



dQ = cud^dfaJdJdiJ 2 - q 2 ) (A.2) 
where cr is a representation dependent constant. These SU(2) color charges satisfy 

dQQ a = 

dQQ a Q? = C 2 5 aP (A.3) 
For fixed Casimir Y2 a Q a Q a = {N 2 — ^)C 2 , we can chose the spherical representation for 

(JAU) 



Q 1 = sin cos 0, Q 2 = sin 9 sin 0, Q 3 = cos9 (A.4) 
for which the measure flA.2j) reads 



dQ = sin OdOdcf) (A.5) 
Equivalently, 
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<? = -^{rr\W)-Y&W)) 

Q 3 = y T n°(^0) (A.6) 

in terms of spherical harmonics. In the spherical representation, we have J dQ = Air, 
E Q Q a Q° = 1 and / dQQ ■ Q = An. 
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